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Abstract 

Unlike the classical exponential relaxation law, the widely prevailing universal 
law with its fractional power-law dependence of susceptibility on frequency cannot 
be explained in the framework of any intuitively simple physical concept. The re- 
sulting constancy of the ratio of the imaginary to the real parts of the complex 
susceptibility, known as the "energy criterion", has a pleasing simplicity but the 
understanding of its origins needs a special theoretical treatment. A fresh light on 
the stochastic nature of the dielectric relaxation has been shed by a novel stochas- 
tic approach introduced in the last decade. Since the theoretical analysis involved 
is rather unfamiliar, the aim of this paper is to give some useful comments and 
suggestions which should help to follow in details the proposed stochastic scheme 
of relaxation leading to the well-known empirical responses. We justify the univer- 
sality of the power-law macroscopic response as well as Jonscher's screening and 
energy criterion ideas, and we give a new basis to the research into the significance 
of relaxation processes. 

Keywords: Dielectric relaxation; screening; energy criterion; heavy-tailed distribu- 
tion. 
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1. Introduction. 



Dielectric relaxation in solids, a process commonly denned as an approach to equilibrium 
of a dipolar system driven out of equilibrium by a step or alternating external electric field, 
represents one of the most intensively researched topics in experimental and theoretical 
physics, see e.g. [1-21]. Wide-ranging experimental information resulting from the studies 
of the basic dielectric characteristics (i.e., the time decay of the depolarization current 
i(t) and the frequency- dependent complex dielectric permittivity e(u) or susceptibility 
x{w) oc e(u) — Cqo, where is the asymptotic value of the dielectric permittivity e{uS) 
at high frequencies) has led to the conclusion that the classical phenomenology of relax- 
ation breaks down in complex materials. It has been found that the Debye behaviour, 
represented by the exponentially decaying relaxation function 

<b(t) = exp(-Upt) (1) 

yielding 



and 



/ \ p / \ debit) , , 

i[t) oc f(t) = ^— = w p exp {—Upt) 

X(w) oc 0» = ( e**f{t)dt = — . (2) 
' p 

(where the constant u> p denotes the loss peak frequency) is hardly ever found in nature and 
that the deviations from it for many dielectrics may be relatively large |], |TU|, [Tj| . A class 
of systems exhibiting the non-Debye relaxation patterns includes various complex ma- 
terials such as supercooled liquids, amorphous semiconductors and insulators, polymers, 
disordered crystals, molecular solid solutions, glasses, etc. 

For a long time effort was being diverted to a purely qualitative representation of the 
shape of the non-Debye dielectric functions <p(t) or </>*(u;) in terms of certain mathematical 
expressions without in any way going into the physical significance of these representa- 
tions. It turns out that all dielectric data are characterised well enough by a few empirical 
functions J| |3], 0, |l^]. The time-domain relaxation data usually are fitted by means of 
the stretched exponential relaxation function 

<j>{t) = exp (-(u p t) a ) with 0<«<1 (3) 

known as the Kohlrausch- Williams- Watts (KWW) function and coincident with the Debye 
case (P for a = 1. The most popular analytical expression applied to the complex 
susceptibility data is given by the Havriliak-Negami (HN) function 

0*M = T, J—, \ — v7 ( 4 ) 

where < a, 7 < 1. For a = 1 and 7 < 1, formula (|]) takes the form known as the Cole- 
Davidson (CD) function; for 7 = 1 and a < 1 it takes the form of the Cole-Cole (CC) 
function, and for a—1 and 7 = 1 one obtains the classical Debye form (El). 
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We repeat that these various model functions are essentially attempts to characterise 
the observed behaviour without in any way indicating the physical mechanisms involved. 
To that extent, we note that the problem of understanding the nature of dielectric relax- 
ation is as yet largely open since we do not have a sufficiently general method of approach 
to the relaxation processes. At the same time we note certain fundamental features of re- 
laxation response which make it clear that we are faced with some very general processes 
into which we require a much deeper insight. 

The HN, CC, CD, and KWW model functions all show the high-frequency limit of the 



response characterised by the fractional power law P, [15 



"x!{uj) = tan ( —J x"( u ) oc uj n ~ x for some < n < 1 (5) 

which has the straightforward consequence of the constant ratio of the imaginary to the 
real components of the dielectric susceptibility x(u;) = x'( u ) ~~ ix"{ u )' 

— j^-r = cot (n—) for uo ^> u v . (6) 
X'(v) V 2 J 

Let us note that this relation implies that the ratio of energy lost per radian to the energy 
stored at the peak is independent of frequency. This very simple and intuitively evident 
property is characteristic of the prevailing universal response of most dielectric materials. 
We call this the "energy criterion" and it is neither known nor predicted by any of the 
accepted theories of relaxation. 

Important for our discussion of relaxation phenomena is the concept of residual loss 
22f which remains after the removal of such strongly frequency-dependent processes as 



the direct-current contribution going as tu~ , or low-frequency dispersion where loss is 
proportional to io n ~ x with n — > 0, or dipolar loss which shows pronounced peaks. The 
residue left is broadly frequency-independent or only weakly dependent over extended 
ranges of frequency which may be approximated by the limit n — > 1. Such a limit 
corresponds to frequency-independent components of the dielectric susceptibility, x'W) 
and x"( a; )' w hh their ratio tending to infinity which is impossible practically but can be 
seen experimentally in various approximations. Figure 1 shows a compilation of data for 
different materials taken from much wider sets for different temperatures. These data are 
representative of a very wide range available in the published literature and are quoted 
here purely by way of examples. The data shown correspond to low or relatively low losses 
so that the real part of the permittivity e'(u) does not vary by a large amount over the 
frequency range in question; and it is therefore possible to normalise the data by dividing 
x"{w) by the corresponding e'(uj), thus defining the loss tangent tan 5 = x"(u;)/e / (u;), 
which gives absolute values of loss for the various samples while retaining the frequency 
dependence of x"{ u )- The essential conclusion, coming from the analysis of the data 
presented in Figure 1, is that the residual loss follows in some cases law (|D while in 
other cases it is essentially "flat" in frequency with minor perturbations. It is interesting 
to note that the addition of 17% of carbon black to polyethylene causes a rise of tan 5 
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by one order of magnitude but hardly any change of the frequency dependence. It is 
also noteworthy that the total range of tan S covers at most three decades for materials 
for which the direct-current conductivity varies by many orders of magnitude. It is not 
profitable to fit these data to any particular law since they are the result of the interplay 
of minor accidental factors superimposed on the prevailing broad "flatness" . A theoretical 
justification for the existence of the fractional power law of type (|5|) with the exponent n 
arbitrarily close to unity would be sufficient for the understanding of all flat losses, any 
minor deviations being purely accidental. 

In our theoretical approach to the analysis of the universal-response characteristics 
(H) we make use of the concept of dipolar screening introduced by Jonscher [23 . In the 



framework of this concept, the effectiveness of screening depends on the relative magni- 
tudes of the thermally-activated-dipole density with an energy W, oc exp(— W/kT), 
and the "critical" density Ni oc kT/fi 2 where [i is the dipole moment of the dipole being 
screened. The theory of dipolar screening predicts that for low dipolar density, Nj, <C Ni, 
the number of dipoles within the field of any one dipole increases almost exponentially, 
so that their behaviour is strongly collective and the etire system behaves in the univer- 
sal manner (|5J). On the other hand, for a high dipolar density, 3> Ni, screening is 
effective, meaning that dipoles do not "see" their neighbours and behave individually in 
a Debye-like manner. This is summarized in Table 1. The situation is in many respects 
similar to the classical screening by charged particles but the onset of screening is much 
more rapid with a fall in the dipolar density N^. 

A considerable effort that has been devoted in the past to finding a theoretical expla- 
nation of the empirically observed results points on the two most widespread and at the 
same time least understood properties of the relaxation responses: 

• the existence of the characteristic property (0) with its limiting form of virtually 
frequency-independent loss; 

• the fact that this property of the relaxation responses is common to a very wide 
range of materials with very different physical and chemical interactions. 

As a consequence, in theoretical attempts to model relaxation it has been commonly 
assumed that the empirical relaxation laws reflect a kind of general behaviour which is 
independent of the details of examined systems. This idea has stimulated the proposal 
of several relaxation mechanisms that differ mainly in the interpretation of the relaxation 
function. In recent attempts to find the origins of the non-Debye relaxation patterns the 
idea of complex systems as the "structures with variations" |24| that are characterised 



through a large diversity of elementary units and strong interactions between them is of 
special importance. The evolution in course of time of physical properties of a complex 



system is nonpredictable or anomalous [25|, and the main feature of all the dynamical 



processes in such a system is their stochastic background. In the framework of statistical 
models (see e.g. [4-6, 8, 9, 11-14, 17-21, 26-31] the fact that the large scale behaviour 
of complex systems shows universality, i.e., that it is to some extent independent of the 
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precise local nature of the considered system, should come as no surprise. Intuitively, one 
expects "averaging principles" like the law of large numbers to be in force. However, it 
turns out to be very hard to make this intuition precise in concrete examples of stochastic 
systems with a large number of locally interacting components. 

The empirical facts stress the need for a completely novel approach to the modelling of 
the dielectric relaxation and, to a certain extent, also mechanical relaxation, photoconduc- 
tion, photoluminescence and chemical reaction kinetics (sharing some common features, 
see The need to understand the connections between the macroscopic property (|5|) 

of the relaxing complex system and the statistical properties of individual molecular or 
dipolar species requires the introduction into the relaxation theory of advanced methods 
of stochastic analysis. As shown by us [8, 32-34], a general formalism of limit theorems 
of probability theory plays an important role in constructing tools to relate the local ran- 
dom characteristics of the complex system to the empirical, deterministic relaxation laws, 
regardless of the specific nature of the system considered. The significance of the present 
paper lies in the fact that no one has clarified in a simple and plausible way, let alone one 
based on a stochastic argument, why the universal relation should exist at all and why 
the residual loss is such slowly variable function of frequency. Our approach provides a 
rigorous explanation for the most widely observed form of frequency dependence of the 
permittivity (or susceptibility) and thereby opens up a new and very powerful way of 
interpreting relaxation phenomena not only in the dielectric context. 

The main objective of this paper is to focus on the approach to relaxation in the 
framework proposed in J33, 34 as consistent with Jonscher's screening and energy crite- 
rion ideas and providing their strict mathematical formulation. We also bring to light the 
spatio-temporal scaling conditions hidden behind all the well-known empirical responses. 
For the reader's convenience, the mathematical details neccessary for the stochastic con- 
struction of the effective representation of a relaxing complex system are followed by 
extended comments. 

In Section 2 we introduce the basic mathematical concept underlying the stochastic 
transition of a complex system from its initially imposed state. We show that the relax- 
ational properties of the entire system can be represented by means of a random effective 
relaxation rate which contains information on the internal stochastic structure of the in- 
vestigated system. In Section 3 we point to the origins of the statistical properties of the 
effective relaxation rate. We discuss the role of limit theorems of probability theory as the 
mathematical technique which allows us to derive the explicit relaxation formulas even 
with rather restricted information on local random properties of the system. In Section 4 
we relate the statistical conditions yielding the well-known empirical responses to the 
spatio-temporal scaling properties of the relaxing complex system, and we show that they 
underlie Jonscher's energy-criterion hypothesis. The last section contains conclusions. 
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2. Transition and survival probabilities 



(i) Let us consider a complex physical system containing identical objects undergoing 
irreversible transitions from state A, imposed at time t = 0, to state B at random instants 
of time. States A and B differ in some physical parameter, so that the transition A — > B 
is defined as the change of this particular parameter (changes in all other parameters may 
also have an influence on the transition). Let us choose one of the objects. Consider the 
conditional probability p(t, dt) that this object will undergo the transition during the time 
interval (t, t + dt) if the transition has not occured before time t, i.e. 

p(t,dt) = Pr(t<9 <t + dt\9>t) (7) 

where 9 is the random waiting time for the transition of the chosen object. The conditional 
probability p(t, dt) defined in (|?p can be expressed in a form 

^ V ' ' Pr(9 > t) 

where Pr(9 > t + dt) and Pr(9 > t) are the survival probabilities, i.e. the probabilities 
that the considered object will remain in state A until time t + dt and t, respectively. The 
survival probability of the object can be expressed as 

Pi(9 >t) = l-p(t) 

where 

p{t) = Pi(9 < t) 

is the waiting time distribution of the object, i.e. the total probability of its transition 
A — > B until time t. One can rewrite the survival probability in the form 

Pr(9 > t) = exp(- J r(s)ds) (8) 
o 



which is dependent on a non-negative quantity r(s) called the intensity of transition [35 
This quantity is time-dependent, in general; and moreover, because of random impacts 
affecting each object, for different objects in the system it can take various values at the 
same instant of time 0. 

The survival probability Pr(6* > t) can be derived if one knows the explicit form of 
the intensity r(s). For the time independent intensity r(s) = bo = const one gets 

Pr(# >t) = exp(-b t) 

what recovers the classical exponential evolutionary law for the object; the value bo de- 
termines the relaxation rate of the transition process. If, on the contrary, the intensity of 
transition is essentially time-dependent, then the evolutionary law for the object is of the 
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nonexponential form. In this case the resulting relaxation rate is not directly given by 
the intensity of transition as it is in the exponential case. In fact, taking into account the 
peculiarity of local random environment of an object, one should assume that, in general, 
the survival probability Pr(# > t) has the form of the weighted average of an exponential 
decay with respect to the probability distribution Fp(b) of the relaxation rate of the object 

oo 

Pr(# > t) = f e~ bt dFp{b). (9) 
o 

In other words, the relaxation rate of the object is the random variable (3 such that the 
total survival probability has the form 

Pr(# >t) = (exp(-/3t)) (10) 

where the mean value (■) is taken with respect to the relaxation-rate probability distribu- 
tion Fp(b). Expressions (Q) and fllCf ) are a generalisation of the Debye relaxation for which 
the relaxation rate (3 takes a constant value, say b , with probability 1, Pr(/3 = bo) = 1. 
Then (3 has a degenerate probability distribution, dFp(b) = 5(b — b )db (where d(-) denotes 
the Dirac delta function), and the integral on the right-hand side of (d) equals 

oo 

/ e- bt 5{b - b )db = exp(-M)- 
o 

In general, in a complex system the probability distribution Fp(b) and, consequently, the 
survival probability Pr(# > t) of an object are of unknown forms. 

From equation (H) the intensity of transition from an initial state, r(t), is related to 
the total survival probability Pr(# > t) of the object as follows 

r(t) = -— lnPr(0 > t). (11) 

Comparing (|TT| ) and @ we obtain the relationship between the time-dependent intensity 
of transition and the relaxation rate distribution 

r(t) = --\nJe- bt dFp(b). 
o 

The above mentioned lack of information about the relaxation rate distribution Fp(b) 
yields that the corresponding intensity r(t) is of unknown form. The explicit form of r(t) 
depends on the characteristics of the random environment around the examined object 
and the rules needed to specify the sets of deterministic and stochastic parameters. 
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(ii) Consider (as before) a system of iV objects, each waiting for transition A — > B for 
some random time. For the ith object, 1 < i < N, let the waiting time be denoted by 9 iN . 
The notation here (i.e. index "iN") emphasizes the impact of the system size iV on the 
behaviour of each individual object that has to be taken into account. The non-negative 
waiting times 9 1N , . . . , 9 NN form a sequence of independent, identically distributed random 
variables. The waiting time distribution Pi^it) = Pr(0 4JV < t) is hence the same for each i 
and is equal to some function F e {t) that may depend on N. 

The objects undergo transition in a certain order that can be reflected in the no- 
tion of order statistics |36j 0(i) < . . . < 9( N ) , which is simply a nondecreasing rearange- 
ment of times 9 1N , . . . , 9 NN . Traditionally, 9^) is called the /-th order statistics of sample 
9 1N , . . . , 9 NN . Note that 9 (i) — min($ 1JV , • • • , 9 NN ) and 0(jv) — max($ 1]V , . . . , 0jvjv). 

For a fixed size iV of the system, the ratio of objects not transformed up to time t is 
equal to 

i-rag. (12) 

where r] N (t) denotes an unknown (random!) number of objects already transformed at 
time t > 0. For I = 0,1, ... N the events {rj N (t) = 1} that up to time t exactly I transitions 
occured in the system can be expressed via order statistics in the following way: 

=0} = {0(i) >t}, 

= {0(0 < t, (z+1) > t} for l = l,...N-l, (13) 
{Vs(t)=N} = {9 {N) <t}, 

and the probability of such events is given by the Bernoulli model 

Pv( VN (t) = I) = ^ [p lN (t)f [l-p lN (t)f- 1 , I = 0,1,... N. (14) 

Thus the random number t] N (t) has the Bernoulli (binomial) distribution B(N,p) with 
parameter p = p iN (t). Consequently, the fraction in ([12]) is a random variable. However, 
it follows from the strong law of large numbers |J7J that for iV large enough 

V*{t) _ / Vif(t) ' 
N ~\ N . 



PlJV (t) = (15) 



i.e. the random nature of the fraction in fll2[) vanishes for large A/". Hence (0) is asymp- 
totically equal to (1 — p iN (t)) which is the survival probability Pr(9 iN > t) of any single 
object: 

1 - ^1 « 1 - p m (t) = Pr(9 lN > t). (16) 

The basic property of the relaxation function <p(t) is its monotonic decrease from 1 
at t = to as t —>■ oo. This is in fact a property of the survival probability and 
therefore (O) has been used as a definition of the relaxation function in several models. 
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Yet (|T2p can be used to describe the evolution of the entire system only if the behaviour 
of the system is represented by any individual object (from those forming the system). 
Unfortunately, this strict condition seems to contradict the idea of complexity of the 
investigated systems, and hence models defining the survival probability of a complex 
system as in (|T6|) do not capture the nature of relaxation phenomena. Moreover, the 
relation flT5| ) holds for any arbitrarily chosen form of the waiting time distribution (e.g. 
F e (t) = for t < 0, = t for < t < 1, and = 1 for t > 1) so that the ratio in flT5|) 
is not uniquely determined. Therefore the attempts leaving the probabilistic analysis of 
the irreversible stochastic transitions in complex systems at this stage can only propose 
the form of distribution Fe{t) fitting the data most exactly, see e.g. [|, ||. They do not 
explain the observed fractional power law (||) indicating strictly limiting properties of the 
survival probability of the initial state of a complex system. 



(iii) The considerations of irreversible stochastic transitions in complex systems show 
that the behaviour of the system as a whole, in general, cannot be attributed to any 
chosen object forming the system |6], || |^, Crucially relevant to this statement is 
the question of a proper mathematical construction of an "averaged" imaginary object 
respresenting the entire system. As it follows from (§), the intensity of transition from 
an initial state for any (real or imaginary) object depends on its survival probability. Let 
us denote the survival probability of an imaginary object representing the system as a 
whole by Pr(6 l jv > t) and by 9n the effective waiting time for the entire system. Since 
the considered imaginary object represents the whole system, the survival probability 
Pr(0jv > t) is the probability that the transition of the system as a whole from its 
initial state (imposed by external constraints at t — 0) has not happened prior to a 
time instant t. This mathematical quantity is defined as the probability that there is no 
transition occuring in the system up to time t so that 

Pr(6W >t) = Pr{rj N {t) = 0). 

By means of the Bernoulli scheme ( ]T4]) with 1 = 0, the survival probability of the system 
is the product of N factors, each asymptotically equal to (0), see equation (|T6|): 

Pr((9jv > t) = Pv( VN (t) = 0) = (l- PlN (t)) N . 

On the other hand, by means of the order statistics {rj^it) = 0} = {9^ > t}, see equation 
(|T3|); and the survival probability of the entire system 

Pr(0;v >t) = Pv( VN (t) = 0) = Pr(0 (1) >t) = Pr(min(0 1JV , . . . , 9 NN ) > t) 

is just the probability that the first passage of the system from the initial state has not 
happened before time t. 

Since the waiting times 9in, . . ■ ,9^n are independent and identically distributed, we 
have 

Pr^ >t) = Pr(9 1N >t)... Pr(9 NN > t). 
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Assuming fllUp for each 6in, i-e. 

Pr(^jv > *) = (eM-M) 

where random variables Pin, ■ ■ ■ , Pnn are also independent and identically distributed, 
we obtain that the survival probability of the system is of the form analogous to (PH): 

Pr(0W>*) = (exp(-^)) (17) 

where 

N 

Pn = Y,&n- (18) 

»=i 

The random variable Pn can be considered as the effective relaxation rate representing the 
entire system; and, in accordance with the rate-theory concept (see e.g. |13|), individual 
relaxation rates Pun's are hence some contributions to the total rate Pn- 

In fact, in any complex dielectric system under a weak external electric field only a part 
of the total number N of dipoles, referred to as "active" dipoles, is directly governed by 
changes of the field. The exact number of such dipoles depends on temperature, density 
of the dipolar species in the system, and interactions between them, and it is usually 
unknown. (In some cases it may take the value of the total number of the dipoles in 
the system.) Therefore, there can be only a part of dipoles contributing to the effective 
relaxation rate, and it is reasonable to modify formula (TIB]) taking some unknown number 
v N instead of exactly N components: 

U N 

Pn = Y,&n- (I 9 ) 

i=i 



3. Origins of the random effective relaxation rate 

The relaxation function <f>(t), which is the probability that the initially imposed state of 
a macroscopic system survives by time t, is well known to be expressed as the weighted 
average of an exponential decay with respect to the distribution of the effective relaxation 
rate Pn 0, |38], ^9|, and its explicit form depends on statistical properties of the rate. 
In view of (|17D , we can use the survival probability Pr(#jv > t) as a definition of the 
relaxation function and we obtain in this case that 

4>{t) = (exp(-P N t)) (20) 

where Pn is given by (|T9|) . 

The evaluation of Un, the number of components in (|H]), which is in agreement with 



the nature of the relaxation phenomenon, has been presented in |33], |34|. It is based on 
the following idea: 
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Depending on the strength of screening, the active dipoles in the complex system may 
"see" each other to some extent. The collective behaviour of active dipoles results in the 
appearance of mesoscopic cooperative regions contributing to the macroscopic relaxation 



process p], The number of such mesoscopic regions in the system is determined 



by their sizes Mi, M 2 , . . . and by the number of the active dipoles in the system [33]. 

The effective relaxation rate (3^ consists of the contributions (3jN of all mesoscopic 
cooperative regions 

h = TA^- ( 21 ) 

3=1 



Similarily, the relaxation rate Pjn of the jth mesoscopic region is the sum of the contri- 
butions PiN of all active dipoles over the region [13, 32-34] and hence it is equal to 

Mi+...+Mj 

= E Ajv- (22) 

i=Afi+...+Mf_i+l 

As mentioned in section 2 (i), the individual relaxation rate, denoted here by reflects 
the random intra-cluster dynamics, i.e. statistical properties of the interactions of the zth 
active dipole with some inactive neighbours forming around it a cluster of size iVj. 



Comparing (]L9|) and (pT|) with fij^ of the form (|22|) one obtains the explicit formula 
for i/jv, the number of the contributions to the total relaxation rate, in the framework of 
the proposed approach. For the sake of this paper it is not neccessary to cite the entire 
complicated formula derived (for details see J33|). It is only important to note that z/jv 
is fully determined by N, the system size; by the cluster sizes Ni,N 2 ,...; and by the 
cooperative- region sizes Mx,M 2 , .... The number of the dipoles directly engaged in the 
relaxation process, their location, interaction range, relaxation rates, and all the quantities 
defined by them are random. Their precise stochastic characteristics determining the 
relaxation response of the system are unknown, unfortunately. However, as a rule, the 
relaxing systems consist of a large number of dipoles so that the distribution of (3n can 
be satisfactorily approximated by the weak limit 

j3 = lim j3 N . (23) 

N— »oo 

In practice, even N ~ 10 6 can suffice to replace adequately (3n in (|20D by the limit (3. 

On the basis of the limit theorems of probability theory the distribution of the effective 
relaxation rate (3 in (p3|) and the resulting form of the relaxation response (|20D in time or 
the corresponding susceptibility in frequency domain 



*V) oc 0* (a,) = J e-^ -^r(t) dt 

o ^ ' 

can be derived even with rather restricted information on properties of the micro/mesoscopic 



statistical levels in the system. It appears [33. 34] that each of the well-known dielectric 
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responses can be obtained in this way if appropriate conditions are imposed on the ba- 
sic random quantities: the cluster sizes iVj, the sizes Mj of the cooperative regions of 
correlated clusters, and the individual active-dipole relaxation rates fax- 



4. Foundations of the universal response 

Table 2 collects the statistical properties shown Q to underlie the KWW, HN, CC, 
CD, and Debye relaxation responses, given by formulas (|3]) and (|4]). As it follows from the 
table, the empirical responses are realized by various stochastic schemes. The conditions 
imposed on Ni, Mj, and /3iN here take one of two forms; either the distribution of the 
quantity considered is of finite mean value or it has a heavy tail. Let us explain that the 
distribution of a nonnegative random variable, say X, has a heavy tail if for some < a < 1 
the tail function Pr(X>a;) ~ x~ a for large x, i.e., 

Pr(X>x) 

lim — i = const > 0. (24) 

x^oo X~ a 



Classical example of the heavy-tailed distribution is the Pareto distribution [HO 



Pr(X>a;) = x>0 

v ; l + (Ax) c 

with the shape parameter < c < 1 and any positive scale constant A. (In this case the 
heavy-tail exponent a = c.) Another distribution possesing this property is a generalisa- 



tion of the Pareto example; namely, the Burr law p0[ 



Pr(X>x) = - t- r-r-r . £>0 

v ; (i + (Ax) c ) d1 

with the shape parameters c, d > such that < cd < 1 and any positive scale constant A. 
(Here the heavy-tail exponent a = cd.) More difficult but interesting heavy-tailed dis- 
tributions are the completely asymmetric Levy-stable laws |36j defined by their Laplace 
transforms C(X; t) = exp(— (At) c ) with < c < 1 and A > 0. (In this case the heavy-tail 
exponent a = c.) 

Condition ( |2"4"D applied to any random variable X expresses the following scaling prop- 
erty of the magnitude represented by X: 

Pr(X > Cx) w C- a Pr(X > x), x -> oo, (25) 

for any fixed constant C > 0. (If there exists the corresponding probability density w(x), 
the condition (|25| ) yields w(Cx) ~ C~ a ~ 1 w(x) for large x.) Hence, in the presented 
approach the heavy-tail property is directly related to the spatial (if referred to iVj and 
Mj, the cluster and the cooperative- region sizes) or temporal (if referred to the relaxation 
rate /%jv) scaling properties of the system. Comparing the non-Debye responses with the 
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Debye one, as done in Table 3, one can see that the difference between them lies in the 
presence of such a scaling in the system. Namely, in the Debye case, that of the individual 
behaviour of the active dipoles, there is no property of this form. By contrast, spatial 
or/and temporal scaling can be found for other responses. Moreover, at least one of the 
two variates Mj and finsr has a heavy-tailed distribution in the non-Debye cases. Since, 



as equation ( [22] ) shows, the correlated-cluster relaxation rate (3j N includes both those 
quantities; its distribution has also a heavy tail, "producing" a hierarchy of mesoscopic 
relaxation rates; i.e. the same proportion of smaller or larger mesoscopic contributions 
to the effective response no matter the scale at which one is looking at the relaxation 
rate distribution. This property follows from both the stochastic properties of the active 
dipoles generated by the intra-cluster dynamics and the strength of screening yielding the 
range of interactions between them. 

As a consequence of the proposed approach, the origins of the macroscopic energy 
criterion coefficient can be pointed out. As we have already explained, the heavy-tailed 
distribution of any (micro/meso/macroscopic) relaxation rate f3 reflects the scaling prop- 
erty of the rate 

Pr(/3 > Cb) w C- a Pr(/3 > b), b -> oo, (26) 

for any fixed constant C > (where a is the heavy-tail exponent). As in flIU| ) and (|TTD, 
we have Pr(6* > t) = (exp (—/%)) so that the properties of the micro/meso/macroscopic 
relaxation rates determine the survival probabilities of the microscopic (active dipole), 
mesoscopic (cooperative region), and macroscopic (entire system) objects, respectively; 
thus, by analogy to (^), defining the function (p{t) = (exp(— fit)), for simplicity referred 
to as the micro/meso/macroscopic relaxation function. 

The asymptotic behavior (|26|) of the relaxation-rate distribution at large b is connected 
with the short-time behavior of the associated relaxation function <f>(t) [[36|]. Namely, at 

the origin t — ► the response function f(t) = — takes the form 

f(t) oc t°- l L{t) (27) 

where L(t) is a function slowly varying at t — (i.e., L(Ct)/L(t) —+ 1 as t — > for 
any fixed positive constant C). Then, it can be shown |^T[ that the short-time property 
(|27| ) of the response function f(t) corresponds to the following asymptotic behavior of 

oo 

x(cj) oc 4>*{oj) = J e~ tujt f(t)dt for the high-frequency region: 
o 

xM = x'H-^"H«M" a i(lM for uj^oo. (28) 
Property ( p!f) leads straightforwardly to the constant ratio 

- cot 7]— tor oj S> uj p 



x'M V 2 
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with 7] = 1 — a. This result, interpreted as the (micro/meso/macroscopic) energy criterion, 
is consistent with the energy criterion (Q) if only 77 is equal to the macroscopic high- 
frequency power-law exponent n. Summing up, the heavy-tail property of the relaxation 
rate with the heavy-tail exponent a leads to the energy criterion with the characteristic 
constant 77 = 1 — a. 

Table 4 contains the detailed discussion of the energy criterion at the micro/meso/ 
macroscopic levels for the cases related to the KWW, HN, CC and CD empirical re- 
laxation responses (^) and (|4j). Observe that in each particular response the parameter 
77 followed from the mesoscopic-relaxation-rate tail is equal to the macroscopic energy- 
criterion coefficient n, and hence, that the high-frequency universal laws (||) and (f|) and 
its parameter n are determined by the heavy-tail properties of the mesoscopic (correlated- 
cluster) relaxation rate distributions. The above proves Jonscher's hypothesis introduced 
in his screening theory of relaxation and stating that in order to observe property (^) at 
the macroscopic level, the same property must characterise the behaviour of the relaxing 
dipoles. As shown in Table 4, in general, the energy-criterion parameter agrees with n 
from equation (H) on the mesoscopic level, only. 

5. Conclusions 

To our knowledge, this is the first attempt to relate rigorously the widely observed uni- 
versal macroscopic laws to the micro/mesoscopic statistical characteristics of the complex 
system, thereby establishing a basis for the understanding of the stochastic origins of the 
relaxation phenomenon. The essential element which was required for the derivation of 
these results is the strict probabilistic formalism in terms of which the dipolar screening 
and energy criterion concepts can be expressed. 

In terms of our analysis, the conditions for the universal relaxation response may be 
stated as follows: 

(i) at the microscopic level, a random number of the active dipoles, those that follow 
changes of the external field, is selected; their individual relaxation rates are de- 
termined by the interactions of the active dipoles with inactive neighbours forming 
random-sized clusters around them; 

(ii) at the mesoscopic level, correlated-cluster regions of sizes depending on the strength 
of screening appear; the collective rate of relaxation of the active dipoles in such 
a mesoscopic cooperative region becomes correlated with the number of the active 
dipoles in the region and with the stochastic properties of their individual relaxation 
rates, see equation fl22f) ; the mesoscopic correlated-cluster relaxation rate combines 
both the spatial and the temporal scaling properties of the system; 

(iii) at the macroscopic level, averaging over the number of effective contributions and 
their rates leads to the universal relaxation for the entire system, giving law ([5]); the 
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sufficient condition of this is that the high-frequency energy criterion is satisfied at 
the mesoscopic level; the mesoscopic energy criterion is provided by the heavy-tail 
property of the correlated-cluster relaxation rates. 
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Table 1. 

Dipolar screening schemes. 



Table 1 



Condition 


Screening 


Physical picture 


Consequences 


N d » TVi 


effective 


Dipoles do not 
"see" one another 


Individual 
behavior 


Debye-like 
response 


N d « jVi 


ineffective 


Dipoles interact 
strongly 


Collective 
behavior 


Universal 
response 
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Table 2. 

Stochastic origins of the empirical relaxation responses. 



Table 2 



Conditions for: 


Relaxation 
response 


Energy- 
criterion 
parameter 


active-dipole 
relaxation 
rate pin 


cluster 
size Ni 


cooperative 
-region 

k 71/7" 

size Mj 


0<n<l 


heavy tail 
with a = a 


heavy tail 
with a = a 


hacnnr Toil 

with = 7 


HN 

0<a< 1 
0<7<1 


n = 1 — 07 


heavy tail 
with a = a 


heavy tail 
with a = a 


(Mj) <oo 


cc 

0<a<l 
7 = 1 


n = 1 — a 


((3 iN )<oo 


(iVi)<oo 


heavy tail 
with = 7 


a = l 

0<7<1 


n = l— 7 


heavy tail 
with a = a 


(iVi)<oo 


{Mj)«x> 


KWW 

0<a< 1 


n = l — a 


{(3 iN )<oo 


(iVi)<oo 


(Mj) <oo 


Debye 

a = l 
7 = 1 


energy 
criterion 

not 
fulfilled 
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Table 3. 

Comparison of the non-Debye relaxation responses to the Debye one. 



Table 3 



Non-Debye responses 

(collective behavior of the dipoles): 


Debye response 

(individual behavior of the dipoles): 


The relaxation rates of the dipoles 
responding to changes of the external 
electric field, 


The relaxation rates of the dipoles 
responding to changes of the external 
electric field, 


or 


and 


the cluster sizes, 


the cluster sizes, 


or 


and 


the sizes of correlated-cluster regions 

have heavy-tailed distributions, 
and, as a result, 
infinite expected values. 


the sizes of correlated-cluster regions 
(if any) 

have distributions 
with finite mean values. 


(Note that the HN response is the only case in which all 
three distributions have heavy tails.) 


(Note that the dipoles do not have to respond with the 
same rate but the rates are of limited distribution.) 


Spatial and/or temporal scaling is present at 
the different statistical levels of the system. 


Neither spatial nor temporal scaling appears 
at any level. 
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Table 4. 

Discussion of the energy criterion at the micro/meso/macroscopic levels. 



Table 4 



Universal 
response 


Relaxation rate 


X"(w)/x'(w) = cot (7777/2) 


PiN 
(micro) 


PjN 

(meso) 


P 

(macro) 


micro 


meso 


macro 


HN 


heavy tail 
with a = a 


heavy tail 
with a = 07 


heavy tail 
with a = «7 


7] = 1 — a ytz n 


7] = 1 — 07 = n 


CC 
KWW 


heavy tail 
with a = a 


heavy tail 
with a = a 


heavy tail 
with a = a 


7] = 1 — a = n 


CD 


((3 iN )<oo 


heavy tail 
with = 7 


heavy tail 
with = 7 


does not 
hold 


7] = 1— 7 = n 
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Figure captions 
Figure 1. 

Typical examples of the frequency dependence of the residual loss, after substraction of 
dc component, low-frequency dispersion and strong dipolar peaks, plotted logarithmi- 
cally against frequency, for a range of dielectric materials. The data cover more than 
ten decades of frequency, and the prevailing trend is a "flat" or frequency-independent 
loss with only minor deviations, for instance the slope —0.2 over six decades for pure 
polyethylene. The tan 8 values cover a range of barely three decades between the lowest 
and highest losses; and the surprising feature is the very slight effect of the addition of 
conducting species (like 12% graphite to polyethylene raising the loss by a factor of 3, 
and 17% of graphite by a further factor of 3) with very little change of the frequency 
dependence. Likewise, the addition of 55% of graphite to boron nitride does not raise 
the extremely low loss of the system. Most of the data quoted here are taken from 
much wider ranges for variable temperature and for ranges of composition. The evident 
trends are the general "flatness" of the losses in frequency which is not compatible with 
any known mechanism and that irrespective of the absolute level of loss. Detailed in- 
formation: Pb2Sco.9Tao.9Tio.1O6 and Pb2ScTaC>6 are ferroelectric ceramics, from A. Isnin 
and A.K. Jonscher, Ferroelectrics, 210 (1998) 47. Pure polyethylene (PE) showing nine 
decades of frequency with a very slowly varying tan<5 in the range 3 x 10~ 4 - 2 x 10~ 3 
with a frequency dependence proportional to /~ ' 2 over most of the middle range; also 
shown are data for polyethylene with admixtures of 12% and 17% of graphite; from 
D.S. McLachlan and M.B. Heaney, Phys. Rev. B 60 (1999) 12746. Single crystal CaTi0 3 
over a range of temperatures with tan<5 = 5 — 9 x 10~ 4 , show extremely small values of 
1-n = 0.000325 to 0.000625; from B.S. Lim, A.V. Vaysleyb, and A.S.Nowick, Applied 
Physics A 56 (1993) 8. Alpha-poly-vinylidene di-fluoride (PVDFa) data from A.K. Jon- 
scher and G. Menegotto, IEEE Trans. Dielectrics and EI 7 (2000) 303. Percolation system 
graphite-boron nitride (45%BN+55%C) from J.J. Wu and D.S. McLachlan, Phys. Rev. B 
56 (1997) 1236; Phys. Rev. B 58 (1998) 14880. Zeolite plus NLO-active para-aniline from 
M. Wuebbenhorst, private communication. 
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Residual Loss, tan 8 
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